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ABSTRACT 

If the Universe contains at least one inflationary domain with a sufficiently large 
and homogeneous scalar field 0, then this domain permanently produces new infla- 
tionary domains of all possible types. We show that under certain conditions this 
process of the self-reproduction of the Universe can be described by a stationary dis- 
tribution of probability, which means that the fraction of the physical volume of the 
Universe in a state with given properties (with given values of fields, with a given 
density of matter, etc.) does not depend on time. This represents a strong deviation 
of inflationary cosmology from the standard Big Bang paradigm. 
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The first models of inflation were based on the standard assumption of the Big Bang theory 
that the Universe was created at a single moment of time in a state with the Planck density, and 
that it was hot and large (much larger than the Planck scale M" 1 ) from the very beginning. The 
success of inflation in solving internal problems of the Big Bang theory apparently removed the 
last doubts concerning the Big Bang cosmology. It remained almost unnoticed that during the 
last ten years the inflationary theory has broken the umbilical cord connecting it with the old Big 
Bang theory, and acquired an independent life of its own. For the practical purposes of description 
of the observable part of our Universe one may still speak about the Big Bang. However, if one 
tries to understand the beginning of the Universe, or its end, or its global structure, then some 
of the notions of the Big Bang theory become inadequate. 

For example, already in the first version of the chaotic inflation scenario |l| there was no 
need to assume that the whole Universe appeared from nothing at a single moment of time 
associated with the Big Bang, that the Universe was hot from the very beginning and that the 
inflaton scalar field <fi which drives inflation originally occupied the minimum of its potential 
energy. Later it was found that if the Universe contains at least one inflationary domain of a 
size of horizon ('/t-region') with a sufficiently large and homogeneous scalar field <j), then this 
domain will permanently produce new /i-regions of a similar type. In other words, instead of 
a single Big Bang producing a one-bubble Universe, we are speaking now about inflationary 
bubbles producing new bubbles, producing new bubbles, ad infinitum. In this sense, inflation 
is not a short intermediate stage of duration ~ 10~ 35 seconds, but a self-regenerating process, 
which occurs in some parts of the Universe even now, and which will continue without end. The 
most striking realization of this scenario occurs in the context of chaotic inflation 0, but the 
basic features of this scenario remain valid in old inflation 0], new inflation ||, |5|] and extended 
inflation as well 0. 

Thus, recent developments of inflationary theory have considerably modified our cosmological 
paradigm |7|. Now we must learn how to formulate physical questions in the new context. For 
example, in a homogeneous part of the Universe there is a simple relation between the density of 
matter and time. However, on a very large scale the Universe becomes extremely inhomogeneous. 
Its density, at the same 'cosmic time', varies anywhere from zero to the Planck density. Therefore 
the question about the density of the Universe at the time 10 10 years may not have any definite 
answer. Instead of addressing such questions we should study the distribution of probability of 
finding a part of the Universe with given properties, and find possible correlations between these 
properties. 

It is extremely complicated to describe an inhomogeneous Universe and to find the corre- 
sponding probability distribution. Fortunately, there exists a particular kind of stationarity of 
the process of the Universe self-reproduction which makes things more regular. Due to the no- 
hair theorem for de Sitter space, the process of production of new inflationary domains occurs 
independently of any processes outside the horizon. This process depends only on the values of 
the fields inside each ^,-region of radius H~ l . Each time a new inflationary /i-region is created 
during the Universe expansion, the physical processes inside this region will depend only on the 
properties of the fields inside it, but not on the 'cosmic time' at which it was created. 
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In addition to this most profound stationarity, there may also exist some simple stationary 
probability distributions which may allow us to say, for example, what the probability is of finding 
a given field at a given point. To examine this possibility one should consider the probability 
distribution P c (<f>,t\x), which describes the probability of finding the field at a given point at a 
time t, under the condition that at the time t — the field at this point was equal to x HI H- 
The same function may also describe the probability that the scalar field which at time t was 
equal to 0, at some earlier time t = was equal to x- 

The probability distribution P c has been studied by many authors, see e.g. 0, @-||. Our 
investigation of this question has shown that in all realistic inflationary models the probabil- 
ity distribution P c {4>, t\x) is n °t stationary [0, M. The reason is very simple. The probability 
distribution P c is in fact the probability distribution per unit volume in comoving coordinates 
(hence the index c in P c ), which do not change during expansion of the Universe. By considering 
this probability distribution we neglect the main source of the self-reproduction of inflationary 
domains, which is the exponential growth of their volume. Therefore, in addition to P c , we in- 
troduced the probability distribution P p (<f),t\x), which describes the probability to find a given 
field configuration in a unit physical volume ||. In the present paper (see also [H| for a more 



detailed presentation) we will show that under certain conditions the stationary probability distri- 
bution P p (4>, t\x) does exist, and a typical relaxation time during which the distribution P p (<p, t\x) 
approaches the stationary regime is extremely small. 

First of all we should remember some details of stochastic approach to inflation. Let us con- 
sider the simplest model of chaotic inflation based on the theory of a scalar field minimally 
coupled to gravity, with the effective potential V(4>). If the classical field is sufficiently homo- 
geneous in some domain of the Universe, then its behavior inside this domain is governed by the 
equations 

+ 3#0 = -dV/d(f) , (1) 

H2 + ? = ik & + ™>) ■ < 2 > 

Here H = a/a, a(t) is the scale factor of the Universe, k = +1, —1, or for a closed, open or flat 
Universe, respectively. M p is the Planck mass, which we will put equal to one in the rest of the 
paper. 

Investigation of these equations has shown that for many potentials V{<p) (e.g., in all power- 
law V(4>) ~ 0™ and exponential V{<p) ~ e°"^ potentials) there exists an intermediate asymptotic 
regime of slow rolling of the field and quasi-exponential expansion (inflation) of the Universe 



0. At this stage the Hubble parameter is H((f)) = y87rV(0)/3. In the theories V(0) ~ 0" 
inflation ends at e ~ n/12. In the theory with V((j)) ~ e a ^ inflation ends only if we bend the 
potential at some point e ; for definiteness we will take e = in this theory. 

Inflation stretches all initial inhomogeneities. Therefore, if the evolution of the Universe were 
governed solely by classical equations of motion, we would end up with an extremely smooth 
Universe with no primordial fluctuations to initiate the growth of galaxies. Fortunately, new 
density perturbations are generated during inflation due to quantum effects. The wavelengths of 
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all vacuum fluctuations of the scalar field grow exponentially in the expanding Universe. When 
the wavelength of any particular fluctuation becomes greater than H~ x , this fluctuation stops 
oscillating, and its amplitude freezes at some nonzero value 5<p(x) because of the large friction 
term 3H(p in the equation of motion of the field 0. The amplitude of this fluctuation then remains 
almost unchanged for a very long time, whereas its wavelength grows exponentially. Therefore, 
the appearance of such a frozen fluctuation is equivalent to the appearance of a classical field 
5<j)(x) that does not vanish after averaging over macroscopic intervals of space and time. 

Because the vacuum contains fluctuations of all wavelengths, inflation leads to the creation 
of more and more perturbations of the classical field with wavelengths greater than H~ l . The 
average amplitude of such perturbations generated during a time interval H~ l (in which the 
Universe expands by a factor of e) is given by 

(3) 

The phases of each wave are random. Therefore, the sum of all waves at a given point fluctuates 
and experiences Brownian jumps in all directions in the space of fields. 

The standard way to describe the stochastic behavior of the inflaton field during the slow- 
rolling stage is to coarse-grain it over /i-regions and consider the effective equation of motion of 
the long- wavelength field ||: 

dt ~ 3H(<j>) 2* U ' ■ { ' 

Here H = ^j8nV/3 , £ (t) is the effective white noise generated by quantum fluctuations, which 
leads to the Brownian motion of the classical field 4>. 

This Langevin equation leads to two stochastic equations for the probability distribution 
P c (4>,t\x)- The first one is called the backward Kolmogorov equation, 

—dT~ = a* [-*r ^ Pc(0 ' tlx) ) - mx)dx~ m ' lx) ■ (5) 

In this equation one considers the value of the field <ft at the time i as a constant, and finds the 
time dependence of the probability that this value was reached during the time t as a result of 
diffusion of the scalar field from different possible initial values x = 0(0). 

The second equation is the adjoint to the first one; it is called the forward Kolmogorov 
equation, or the Fokker-Planck equation, 

dPMAx) 1 d (HWy) d fH^l x V'((f>) A 

at = 2dj \—2^T^—2^ mAx) ) + mW) mAx) ) ' (6) 

This equation was derived in M, see also B. 
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One may try to find a stationary solution of equations (|6|), assuming that 9Pc ^ x ^ = 0. 
The simplest stationary solution (subexponential factors being omitted) would be 

This function is extremely interesting. Indeed, the first term in ([/]) is equal to the square of the 



Hartle-Hawking wave function of the Universe JlTJ, whereas the second one gives the square of 
the tunneling wave function [I2| ! 



At first glance, this result gives a direct confirmation and a simple physical interpretation 
of both the Hartle-Hawking wave function of the Universe and the tunneling wave function. 
However, in all realistic cosmological theories, in which V{4>) = at its minimum, the Hartle- 
Hawking distribution exp ( 8V(<j>) ) * s no ^ norma hzable. The source of this difficulty can be easily 
understood: any stationary distribution may exist only due to compensation of the classical flow 
of the field downwards to the minimum of V(0) by the diffusion motion upwards. However, 
diffusion of the field discussed above exists only during inflation. Thus, there is no diffusion 
motion upwards from the region < <p e . Therefore all solutions of equation (|B|) with the proper 
boundary conditions at = e (i.e. at the end of inflation) are non-stationary (decaying) [[J. 

The situation with the probability distribution P p is much more interesting and complicated. 
As was shown in its behavior depends strongly on initial conditions. If the distribution P p was 
initially concentrated at < 0*, where 0* is some critical value of the field, then it moves towards 
small in the same way as P c , i.e. it cannot become stationary. On the other hand, if the initial 
value of the field is larger than 0*, the distribution moves towards larger and larger values of 
the field 0*, until it reaches the field P , at which the effective potential of the field becomes of 
the order of Planck density M p (we will assume M p = 1 hereafter), where the standard methods 
of quantum field theory in a curved classical space are no longer valid. 



Some further steps towards the solution of this problem were made by Nambu and Sasaki [13] 



and Mijic fll3 |. Their papers contain many important results and insights. However, Mijic |E3| 
did not have a purpose to obtain a complete expression for the stationary distribution P p (4>, t\x)- 
The corresponding expressions were obtained for various types of potentials V(4>) in [|13[]. Unfor- 
tunately, according to fl3|| , the stationary distribution P p ((p,t\x) is almost entirely concentrated 
at ^> P , i.e. at V(<f>) ^> 1, where the methods used in |13| are inapplicable. 

We will continue this investigation by writing the system of stochastic equations for P p . These 
equations can be obtained from eqs. (|5|), @ by adding the term 3HP P , which appears due to the 
growth of physical volume of the Universe by the factor 1 + 3H(<f)) dt during each time interval 
dt 0-1T1, 0: 



dP p _ i d ( Hnt) d ^(0) x H0) p ) , 3HU)P (Q) 

~dt ~ 290 \~^r~ 90 y~^r Pp ) + swm Pp ) + 3H{(l))Pp ■ (9) 
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To find solutions of these equations one must specify boundary conditions. The boundary condi- 
tions at the end of inflation follow from the conservation of the probability flux |H| : 



d_ 

30 



3 V^_ 

4 V 



P P (<Pe) , 



d_ 

dx 



( 



x=4> e 



. 



(10) 



The last boundary condition is especially interesting, since it indicates that at least in the vicinity 
of x = 0e the probability distribution with respect to x looks as a square of the tunneling 
wave function. However, we have found solutions of (H), @ to be rather stable with respect 
to modification of these boundary conditions, whereas the conditions at the Planck boundary 
(f) = (f) p do play a very important role. 

Investigation of these conditions poses many difficult problems. First of all, our diffusion 
equations are based on the semiclassical approach to quantum gravity, which breaks down at 
super-Planckian densities. Secondly, the shape of the effective potential may be strongly modified 
by quantum effects at ~ <p p . Finally, our standard interpretation of the probability distribution 
P c and P p breaks down at V(4>) > 1, since at super-Planckian densities the notion of a classical 
scalar field in classical space-time does not make much sense. 

However, these problems by themselves suggest a possible answer. Inflation happens only in 
theories with very flat effective potentials. At the Planck density nothing can protect the effective 
potential from becoming steep. Hence, one may expect that inflation ceases to exist at <ft > §v> 
which leads to the boundary condition 



There is also another, much more general reason to expect that inflation kills itself as the 
potential energy density approaches the Planck density V ~ 1. Indeed, the amplitude of fluc- 
tuations of the scalar field generated during the typical time 5t = H~ x is given by H/2ir, and 
their typical wavelength at that time is 0{H~ 1 ). This means that the energy density associated 
with the gradients of these perturbations is of the order of if 4 ~ V 2 . Thus, in the domains with 
V > 1 the gradient energy density ~ V 2 becomes larger than the potential energy density V. 
This violates one of the basic assumptions necessary for inflation in domains with V > 1. One 
may expect also that large gradients of energy on the scale comparable to the scale of the horizon 
H^ 1 lead to creation of black holes rather than to the permanent self-reproduction of inflationary 
^i-regions. 

Of course, one may argue that all our considerations do not make sense at densities larger 
than the Planck density. When the energy density in any /i-region approaches the Planck density, 
it may no longer be described in terms of classical space-time and should be just thrown away 
from our consideration. In particular, its volume should not be considered as contributing to the 
total volume of the Universe. Thus, such domains should be neglected in our definition of P p . In 
this case the distribution P p for the field will stop moving towards higher values of (ft and will 
approach a stationary regime when this distribution will be shifted towards <fi ~ <p p . 




(11) 



where V(<f> p ) = V( Xp ) = 0(1). 
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What we are saying is even stronger. Even if one makes an attempt to consider the domains 
with V > 1 as a part of classical space-time, many parts of these domains drop out from the 
process of inflation. This means that the total volume of inflationary /i-regions cannot grow as 
fast as e 3Ht . 

We do not know which of these arguments, if any, will survive in the future theory of all 
fundamental interactions. However, all these arguments point out in the same direction: For a 
phenomenological description of stochastic processes in classical space-time one should impose 



boundary conditions of the type of (11) which do not permit penetration of inflation deep into 



the realm of super-Planckian densities. As we will show in |I0| , the exact form of these boundary 



conditions is not very important; most of them allow the same class of solutions as the boundary 



condition (11). Therefore in this paper we will use these boundary conditions, assuming for 



defmiteness that they are imposed at V{(j) p ) = V(x P ) = 1. 

One may try to obtain solutions of equations (|8]), (|9|) in the form of the following series of 
biorthonormal system of eigenfunctions of the pair of adjoint linear operators (defined by the left 
hand sides of the equations below): 

oo 

^(<Mx) = E eAst ^(xK(<« • C 12 ) 



s=l 



Indeed, this gives us a solution of eq. (§) if 
and 



The orthonormality condition reads 



ips{x) Kj{x) dx = 5 sj . (15) 



In our case (with regular boundary conditions) one can easily show that the spectrum of Xj is 
discrete and bounded from above. Therefore the asymptotic solution for P p (<j>,t\x) (i n the limit 
t — > oo) is given by 



P p (0,t| X )=e Al Vi(xW0) .(l + O 



-(Ai-A 2 )i 



(16) 



Here ipi( X ) is the only positive eigenfunction of eq. (|13D, Ai is the corresponding (real) eigenvalue, 
and 7Ti(0) is the eigenfunction of the conjugate operator (0) with the same eigenvalue \±. Note, 
that Ai is the highest eigenvalue, Re (Ai — A2) > 0. This is the reason why the asymptotic 
equation ( |16"D is valid at large t. We have found [10| that in realistic theories of inflation a typical 
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time of relaxing to the asymptotic regime, At ~ (Aj — A2) _1 , is extremely small. Typically it is 
only about a few thousands Planck times, i.e. about 10 -40 sec. This means that the normalized 
distribution 

P p (<M| X )=e- Alt P p (<M|x) (17) 
rapidly converges to the time-independent normalized distribution 

P P (<P\X) = P P (<P,t^ oo\ X ) = MX) 7Tl(0) . (18) 

It is this stationary distribution that we were looking for. Because the growing factor e~ Al * is the 
same for all (and x)i one can nse Pp instead of P p for calculation of all relative probabilities. 
In particular, P p (0|x) gives us the fraction of the volume of the Universe occupied by the field 
0, under the condition that the corresponding part of the Universe at some time in the past 
contained the field x- The remaining problem is to find the functions ipi(x) and 7Ti(0), and to 
check that all assumptions about the boundary conditions which we made on the way to eq. flOp 
are actually satisfied. 



We have solved this problem for chaotic inflation in a wide class of theories including the 
theories with polynomial and exponential effective potentials V(<p) and found the corresponding 
stationary distributions Here we will present some of our results for the theories ^0 4 and 

V e a<i >. 



Solution of equations (|T3|) and ( HP for ipi(x) an d fl"i(0) in the theory |0 4 shows that these 
functions are extremely small at ~ e and x ~ Xe- They grow at large and Xi then rapidly 
decrease, and vanish at = x = 4>p- With a decrease of A the solutions become more and 
more sharply peaked near the Planck boundary. (The functions ipi and 7Ti for the exponential 
potential have a similar behavior, but they are less sharply peaked near P .) A detailed discussion 



of these solutions will be contained in [10J. The eigenvalues Ai corresponding to different coupling 
constants A are given by the following table: 



A 


1 


KT 1 


1(T 2 


1(T 3 


1(T 4 


10~ 5 


lO" 6 


Ai 


2.813 


4.418 


5.543 


6.405 


7.057 


7.538 


7.885 



One can find also the second eigenvalue A2. For example, for A = 10 4 one gets A2 = 6.789. This 
means that for A = 10~ 4 the time of relaxation to the stationary distribution is At ~ (Ai — A2) _1 ~ 
AM' 1 ~ 10" 42 seconds — a very short time indeed. 

Note that the parameter Ai shows the speed of exponential expansion of the volume filled by a 
given field 0. This speed does not depend on the field 0, and has the same order of magnitude as the 
speed of expansion at the Planck density. Indeed, Ai should be compared to 3H(<f>) = 2y67rV(0), 

which is equal to at the Planck density. It can be shown |K| that in the limit A — >• the 

eigenvalue Ai also becomes equal to 2v / 67r ~ 8.681. The meaning of this result is very simple: in 
the limit A — > our solution becomes completely concentrated near the Planck boundary, and Ai 
becomes equal to 3if (0 P ). 

At first glance, independence of the speed of expansion of volume e Alt on the value of the field 



8 



may seem counterintuitive. The meaning of this result is that the domain filled with the field 
gives the largest contribution to the growing volume of the Universe if it first diffuses towards 
the Planckian densities, spends there as long time as possible expanding with nearly Planckian 
speed, and then diffuses back to its original value 0. 

But what about the field which is already at the Planck boundary? Why do the correspond- 
ing domains not grow exactly with the Planckian Hubble constant H((j) p ) = 2\Z&n/3? It happens 
partially due to diffusion and slow rolling of the field towards smaller 0. However, the leading 
effect is the destructive diffusion towards the space-time foam with > P . One may visualize 
this process by painting white all domains with V((fi) < 1, and by painting black domains filled 
by space-time foam with V(<p) > 1. Then each time H~ 1 ((p p ) the volume of white domains with 
~ <pp grows approximately e 3 times, but some 'black holes' appear in these domains, and, as a 
result, the total volume of white domains increases only e 3 <W 2 v^ times. This suggests (by anal- 
ogy with [|T7|) calling the factor df = 3Xi/2\/6tt 'the fractal dimension of classical space-time', 
or 'the fractal dimension of the inflationary Universe'. (Note that df < 3 for A ^ 0; for example, 
df = 2.6 for A = 10~ 5 .) However, one should keep in mind that the fractal structure of the 
inflationary Universe in the chaotic inflation scenario in general is more complicated than in the 
new or old inflation and cannot be completely specified just by one fractal dimension . 



The distribution P p ((f>\x) = , 0i(x) 7r i(0) which we have obtained does not depend on time t. 
However, in general relativity one may use many different time parametrizations, and the same 
physics can be described differently in different 'times'. One of the most natural choices of time 
in the context of stochastic approach to inflation is the time r = In ^fei = / H((p(x,t),t) dt 
||, U . Here a (x, t) is a local value of the scale factor in the inflationary Universe. By using this 
time variable, we were able to obtain not only numerical solutions to the stochastic equations, 
but also simple asymptotic expressions describing these solutions. For example, for the theory 
|0 4 both the eigenvalue Ai and the 'fractal dimension' df (which in this case refers both to the 
Planck boundary at P and to the end of inflation at e ) are given byd/ = Ai~3 — 1.1 y/X, and 
the stationary distribution is 



3 11 

P P (M X ) ~ ex p(~^)) ( v{x) + o,4 - o) • **("* (3 - Al 



1,2 



Note that the first factor coincides with the square of the tunneling wave function [Q! This 
expression is valid in the whole interval from e to P and it correctly describes asymptotic 
behavior of P p (0, r\x) both at x ~ Xe and at x ~ Xp- 

A similar investigation can be carried out for the theory U(0) = V a e a( ^. The corresponding 
solution is 

PMAx) - «p(-s4)) " ' (4) " v ~ y ^ ■ (20) 

This expression gives a rather good approximation for P p {4>, r\x) for all and x- 
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The main result of our work is that under certain conditions the properties of our Universe 
can be described by a time-independent probability distribution, which we have found for theories 
with polynomial and exponential effective potentials. A lot of work still has to be done to verify 
this conclusion, see fL0"| . However, once this result is taken seriously, one should consider its 



interpretation and rather unusual implications. 

When making cosmological observations, we study our part of the Universe and find that in 
this part inflation ended about t e ~ 10 10 years ago. The standard assumption of the first models 
of inflation was that the total duration of the inflationary stage was At ~ 10~ 35 seconds. Thus 
one could come to an obvious conclusion that our part of the Universe was created in the Big 
Bang, at the time t e + At ~ 10 10 years ago. However, in our scenario the answer is quite different. 

Let us consider an inflationary domain which gave rise to the process of self-reproduction 
of new inflationary domains. For illustrative purposes, one can visualize self-reproduction of 
inflationary domains as a branching process, which gives a qualitatively correct description of 
the actual physical process we consider. During this process, the first inflationary domain of 
initial radius ~ within the time splits into e 3 ~ 20 independent inflationary 

domains of similar size. Each of them contains a slightly different field 0, modified both by 
classical motion down to the minimum of V(0) and by long-wavelength quantum fluctuations of 
amplitude ~ H/2ir. After the next time step which will be slightly different for each of 

these domains, they split again, etc. The whole process now looks like a branching tree growing 
from the first (root) domain. The radius of each branch is given by H~ x ] the total volume of 
all domains at any given time t corresponds to the 'cross-section' of all branches of the tree at 
that time, and is proportional to the number of branches. This volume rapidly grows, but when 
calculating it, one should take into account that those branches, in which the field becomes larger 
than (j) p , die and fall down from the tree, and each branch in which the field becomes smaller than 
e , ends on an apple (a part of the Universe where inflation ended and life became possible). 

One of our results is that even after we discard at each given moment the dead branches 
and the branches ended with apples, the total volume of live (inflationary) domains will continue 
growing exponentially, as e Al *. What is even more interesting, we have found that very soon 
the portion of branches with given properties (with given values of scalar fields, etc.) becomes 
time-independent. Thus, by observing any finite part of a tree at any given time t one cannot 
tell how old the tree is. 



To give a most dramatic representation of our conclusions, let us see where most of the apples 
grow. This can be done simply by integrating e Xlt from t = to t = T and taking the limit as 
T — > oo. The result obviously diverges at large T as A^f 1 e AlT , which means that most apples 
grow at an indefinitely large distance from the root. In other words, if we ask what is the total 
duration of inflation which produced a typical apple, the answer is that it is indefinitely large. 

This conclusion may seem very strange. Indeed, if one takes a typical point in the root domain, 
one can show that inflation at this point ends within a finite time At ~ 10~ 35 seconds. This is 
a correct (though model-dependent) result which can be confirmed by stochastic methods, using 
the distribution P c ((j),At\x) [0- How could it happen that the duration of inflation was any 
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longer than 10 seconds? 

The answer is related to the choice between P c and P p , or between roots and fruits. Typical 
points in the root domain drop out from the process of inflation within 10~ 35 seconds. The number 
of those points which drop out from inflation at a much later stage is exponentially suppressed, 
but they produce the main part of the total volume of the Universe. Note that the length of each 



particular branch continued back in time may well be finite |T8] . However, there is no upper limit 
to the length of each branch, and, as we have seen, the longest branches produce almost all parts 
of the Universe with properties similar to the properties of the part where we live now. Since 
by local observations we can tell nothing about our distance in time from the root domain, our 
probabilistic arguments suggest that the root domain is, perhaps, indefinitely far away from us. 
Moreover, nothing in our part of the Universe depends on the distance from the root domain, 
and, consequently, on the distance from the Big Bang. 

Thus, inflation solves many problems of the Big Bang theory and ensures that this theory 
provides an excellent description of the local structure of the Universe. However, after making 
all kinds of improvements of this theory, we are now winding up with a model of a stationary 
Universe, in which the notion of the Big Bang loses its dominant position, being removed to the 
indefinite past. 

The stochastic approach to inflation used in our work has an intermediate position between 
purely classical and purely quantum mechanical approaches to cosmology. In particular, station- 
arity of our probability distribution is closely reminiscent of the time-independence of the wave 
function of the Universe in the Wheeler-DeWitt equation [l9j . We hope that stochastic methods 



may show us a way towards a complete quantum mechanical description of the stationary state 
of the Universe. 

The authors are grateful to D. Linde for the help with computer calculations and to M. 
Mijic, S. Molchanov, A. Starobinsky and A. Vilenkin for valuable discussions. This research was 
supported in part by the National Science Foundation grant PHY-8612280. 
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